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The transverse-field XY model in one dimension is a well-known spin model for which the ground 
state properties and excitation spectrum are known exactly. The model has an interesting phase 
diagram describing quantum phase transitions (QPTs) belonging to two different universality classes. 
These are the transverse-field Ising model and the XX model universality classes with both the 
models being special cases of the transverse-field XY model. In recent years, quantities related 
to quantum information theoretic measures like entanglement, quantum discord (QD) and fidelity 
have been shown to provide signatures of QPTs. Another interesting issue is that of decoherence to 
which a quantum system is subjected due to its interaction, represented by a quantum channel, with 
an environment. In this paper, we determine the dynamics of different types of correlations present 
in a quantum system, namely, the mutual information I{pab), the classical correlations C (pab) 
and the quantum correlations Q (pab), as measured by the quantum discord, in a two-qubit state. 
The density matrix of this state is given by the nearest-neighbour reduced density matrix obtained 
from the ground state of the transverse-field XY model in Id. We assume Markovian dynamics 
for the time-evolution due to system-environment interactions. The quantum channels considered 
include the bit-flip, bit-phase-flip and phase-flip channels. Two different types of dynamics are 
identified for the channels in one of which the quantum correlations are greater in magnitude than 
the classical correlations in a finite time interval. The origins of the different types of dynamics are 
further explained. For the different channels, appropriate quantities associated with the dynamics 
of the correlations are identified which provide signatures of QPTs. We also report results for 
further-neighbour two-qubit states and finite temperatures. 
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I. INTRODUCTION 

The fully anisotropic transverse-field XY model in one 
dimension (Id) describes an interacting spin system for 
which many exact results on ground and excited state 
properties including spin correlations are known [H-Q. 
The Hamiltonian describing the model is given by 
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(1) 



where L denotes the number of sites in the Id lattice, 
a^'^'^ are the Pauli spin operators defined at the lat- 
tice site i, 7 is the degree of anisotropy (—1 < 7 < 1) 
and A the inverse of the strength of the transverse mag- 
netic field in the z direction (A > 0). The Hamiltonian 
is translation invariant and satisfies periodic boundary 
conditions. Two specific cases of the XY model are the 
transverse field Ising model corresponding to 7 = ±1 and 
the isotropic XX model (7 = 0) in the transverse mag- 
netic field. For all values of the anisotropy parameter, 
the Hamiltonian can be diagonalized exactly in the ther- 
modynamic limit L ^ 00 [ll-Q through the successive 
applications of the Jordan- Wigner and Bogoliubov trans- 
formations. For non-zero values of 7, the model exhibits 
a second-order quantum phase transition (QPT) at the 
critical point Ac = 1 separating a ferromagnetic ordered 
phase from a quantum paramagnetic phase. The order 
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parameter for the transition is (cr^), the magnetization in 
the X direction, which has a non-zero expectation value 
for A > 1. The magnetization in the z direction, (cr^), 
is non-zero for all values of A with its first derivative ex- 
hibiting a singularity at the critical point Ac = 1. In the 
interval < I7I < 1, the critical point transition at Ac = 1 
belongs to the Ising universality class. For A 6 [0,1], 
there is another QPT, termed the anisotropy transition, 
at the critical point 7 = 0. The transition belongs to a 
different universality class and separates two ferromag- 
netic phases with orderings in the x and y directions re- 
spectively [IHll- Figured] shows the phase diagram of the 
transverse-field XY chain in the (A, 7) parameter space. 

The transverse-field XY model and its special case, 
the transverse-field Ising model (TIM), are well-known 
statistical mechanical models which illustrate QPTs. A 
QPT occurs at T = and is brought about by tun- 
ing a non-thermal parameter, e.g., exchange interaction 
strength, external magnetic field etc. In recent years, 
quantum information-related measures like entanglement 
, discord [13, [ill and fidelity [T^ have been used as 
indicators of QPTs. A first-order QPT is characterized 
by a discontinuity in the first derivative of the ground 
state energy with respect to the tuning parameter. Simi- 
larly in a second-order QPT, termed a critical point QPT, 
the second derivative of the ground state energy becomes 
discontinuous or divergent at the critical point. Entangle- 
ment and quantum discord (QD) provide different mea- 
sures of quantum correlations in the mixed state of an 
interacting many-body system. The QD has so far been 
defined in the case of a bipartite quantum system. At 
a first-order QPT point, appropriate entanglement mea- 
sures and the QD are known to become discontinuous 
[13, [13 whereas a critical-point QPT is signaled by the 
discontinuity or divergence of the first derivative of the 
quantum correlation measures with respect to the tun- 



2 



^0 



Critical Ising 



■ Ising (y=1) 



Critical XX 



■ XX (7=0) 



Critical Ising 



■ Ising (y=-1) 



Figure 1. (Color online) Criticality of the transverse- field XY 
model in the (A, 7) parameter space. The transverse-field Ising 
model (7 = ±1) has a QCP at Ac = 1. The line A = 1 
represents criticality of the XY model with universality class 
the same as that of the transverse-field Ising model. The 7 = 
line {XX model) represents the anisotropy transition line for 
A G [0, 1] with the transition belonging to a new universality 
class. 



time. The sudden changes in the decay rates of correla- 
tions and their constancy in certain time intervals have 
been demonstrated in actual experiments (26l . \2l\ . 

In this paper, we consider a two-qubit system each 
qubit of which interact with an independent reservoir. 
The density matrix of the two-qubit system is described 
by the reduced density matrix derived from the ground 
state of the transverse-field XY model in Id. A simi- 
lar study has earlier been carried out for the TIM in Id 
[23 | . We investigate the dynamics of the QD as well as the 
classical correlations under Markovian time evolution and 
identify new features close to the quantum critical points 
of the model. The study is further extended to the cases 
of further-neighbour two-qubit states and finite tempera- 
tures. In Sec. ini the calculational scheme for computing 
C (pab) and Q [pab) is described. Wc further discuss the 
methodology for studying the dynamics of the classical 
and quantum correlations for different quantum channels 
representing the system-environment interactions. Sec. 
mil presents the major results obtained and the discus- 
sions thereof. Sec. IIVI contains some concluding remarks. 



ing parameter [15|. The correlations between the differ- 
ent constituents of an interacting quantum system have 
two distinct components: classical and quantum. While 
entanglement provides the most well-known example of 
quantum correlations, the QD captures quantum corre- 
lations which are more general than those measured by 
entanglement [l^, In fact, there are separable mixed 
states which by definition are unentangled but have non- 
zero QD. 

Quantum systems, in general, are open systems be- 
cause of the inevitable interaction of a system with its en- 
vironment. This results in decoherence accompanied by a 
decay of the quantum correlations in composite systems. 
A number of recent studies [T^-H^l explore the dynamics 
of entanglement and QD when the system-environment 
interactions are taken into account. An interesting result 
yielded by such studies is that the QD is more robust than 
entanglement in the case of Markovian (memoryless) dy- 
namics. The QD decays in time vanishing only asymp- 
totically [13, [22|-[25l w hereas entanglement undergoes a 
'sudden death'jlSl Il9|. i.e., a complete disappearance at 
a finite time. The classical correlations, C {pab), and the 
QD, Q {pab), have mostly been computed only for two- 
qubit systems described by the density matrix pab with 
A and B denoting the subsystems. Three general types 
of Markovian time evolution have been observed [2l[: (i) 
C {pab) is constant as a function of time and Q {pab) 
decays monotonically, (ii) C {pab) decays monotonically 
over time till a parametrized time Psc is reached and re- 
mains constant thereafter. Q {pab) has an abrupt change 
in the decay rate at psc, and has magnitude greater than 
that of C {pab) in a parametrized time interval, and (in) 
both C {pab) and Q {pab) decay monotonically. Mazzola 
et al. [22| have obtained a significant result that under 
Markovian dynamics and for a class of initial states the 
QD remains constant in a finite time interval Q <t <i. In 
this time interval, the classical correlations C {pab) decay 
monotonically. Beyond t — t, C {pab) becomes constant 
whereas the QD decreases monotonically as a function of 



II. CLASSICAL AND QUANTUM 
CORRELATIONS: DEFINITION AND 
DYNAMICS 

In classical information theory, the total correlation be- 
tween two random variables A and B is quantified by 
their mutual information I{A,B) = H{A) + H{B) - 
H{A,B). The random variables A and B take on the 
values 'a' and respectively with probabilities given 
by the sets {pa} and {pt}. H{A) = -Ea-PalogaPa and 
H{B) = —^i,Pb^og2Pb arc the Shannon entropies for 
the variables A and B. H{A, B) ~ — },Pa,b log2 Pa,b is 
the joint Shannon entropy for the variables A and B with 
Pafi being the joint probability that the variables A and 
B have the respective values a and h. Generalization to 
the quantum case is straightforward with the density ma- 
trix p replacing the classical probability distribution and 
the von Neumann entropy S{p) ~ —Tt {plog2 p) replac- 
ing the Shannon entropy. The expression for the quantum 
mutual information is given by 

I {pab) ^ S {pa) + S {pb) - S {pab) (2) 

An equivalent expression for the classical mutual infor- 
mation is J{A,B) = H{A) - H{A\B) where the con- 
ditional entropy H{A\B) is a measure of our ignorance 
about the state of A when that of B is known. The equiv- 
alence with the earlier expression for the classical mu- 
tual information can be shown via the Bayes' rule. The 
quantum versions of J{A, B) and I{A, B) arc not, how- 
ever, equivalent because the magnitude of the quantum 
conditional entropy is determined by the type of mea- 
surement performed on the subsystem B to gain knowl- 
edge of its state. Different measurement choices yield 
different values for the conditional entropy. We consider 
von Neumann-type measurements on B represented by 
a complete set of orthogonal projectors, {H^}, corre- 
sponding to the set of possible outcomes k. The final 
state of the composite quantum system after measure- 
ment on B is given by = (l (g) Ilf ) pab {I <8) Ilf ) /pk 
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Figure 2. (Color online) Bit-phase-flip channel: Decay of 
mutual information I {pab) (solid line), classical correlations 
C (pab) (dashed line) and the QD Q (pab) (dot-dashed line) 
as a function of the parametrized time p — 1 — e~^* for 
A = 0.7,7 = 0.7. Also, psc = 0.114. (inset) (i)The variations 
of |ci(p)| (solid line) and |c2(p)| (dashed line) as functions of 
p. The crossing of |ci(p)| with |c2(p)| sets the parametrized 
time point psc- (ii) Decay of mutual information I (pab) (solid 
line), classical correlations C {pab) (dashed line) and quan- 
tum discord Q {pab) (dot-dashed line) as a function of the 
parametrized time p = 1 — e^"' for A — 0.7, 7 — —0.7. 
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where A and B correspond to the two individual qubits 
and z and / are real numbers. The eigenvalues of pab 
are ^ 

Ao = J |(1+C3)+ \/4c2 + (ci -C2)^ 

Ai = i|(l + C3)-\/4cl + (ci-C2)^ 

1 

C2) 



A2 = 4(l 



C3 + Cl + I 



(1 - ca - ci - C2) 



(8) 



with 



Cl 2z + 2/ 
C2 = 2z - 2/ 
C3 = a + d-2b 

Ci = a — d (9) 
The mutual information / {pab) can be written as [T^] 



with pk = Tr ((/ ® nf ) (/ ® nf )) . / denotes the 
identity operator for the subsystem A and pk is the prob- 
ability for obtaining the outcome k. The quantum con- 
ditional entropy is given by 



S {pAB\{^f]) =Y.PkS [pk) 
k 



(3) 



leading to the following quantum extension of the classi- 
cal mutual information 



J {PAB,{IV'^]) = S{pa) - S {pAB\{li^] \ 



(4) 



The projective measurements on the subsystem B re- 
moves all the non-classical correlations between the sub- 
systems A and B. Henderson and Vedral [l^l have quan- 
tified the total classical correlations C {pab) by maximiz- 
ing J (pas, {nf}) w.r.t. {nf},i.e., 



C{pab) 



max (J (pab, {nf})) 



(5) 



The difference between the total correlations / {pab) (Eq. 
^) and the classical correlations C {pab) (Eq. (O) de- 
fines the QD, Q {pab), 



Q (pab) = I (pab) - C {pab) 



(6) 



Though the concept of the QD is well-established, its 
computation is restricted to mostly two-qubit states. For 
the two-qubit X states, analytic expressions for the QD 
can be derived when the density matrix, defined in the 
computational basis {|11), |10), |01), |00)}, has the foUow- 



/ {pab) = S {pa) + S {pb) + ^ A„ log2 A„ (10) 



where 



S (pa) = S {pb) = - 
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1 — C4 



(11) 



With expressions for I {pab) and C {pab) given in Eqs. 
(dni), dn]) and dll), the QD, Q{pab), (Eq. ^) can in 
principle be computed. The difficulty lies in carrying out 
the maximization procedure needed for the computation 
of C {pab)- It is possible to do so analytically when pab 
is of the form given in ( f7)) resulting in the following ex- 
pressions for the QD |28j : 



where 

Qi 

and 



Q [pab) = min {Qi, (92} 



S {pb) - S {pab) - alog2 



(12) 



-61og2 



- d log2 ■ 



A+ log2 A^ 



(13) 



Q2 = S {pb) - S {pab) 

-A_ log2 A_ (14) 

with A± i (1 ± F) and F^ = (a - d)^ + 4 {\z\ + \f\)^ 

For spin Hamiltonians with certain symmetries, the 
two-spin reduced density matrix pij has a form similar 
to that given in Eq. ([7]) (the qubits A and B now repre- 
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sent the spins at sites i and j respectively) where the ma- 
trix elements of pij can be written in terms of single-site 
magnetization and two-site spin correlation functions. In 
the case of the XY model in a transverse-field, the ma- 
trix elements of the two-site reduced density matrix are 

iiini 



d 
b 

Z : 
/ 



2 



4 



(1 - K-D) 



(15) 



The finite temperature single-site magnetization (cr^) in 
the case of the transverse-field XY model is given by 

aiiii 



) tanli 



1 + A cos ( 



2 / 



where 



(7A sin (f))" + (1 + A cos ( 



(16) 



(17) 



describes the energy spectrum and P = with k being 
the Boltzman constant and T the absolute temperature. 
The spin-spin correlation functions are obtained from the 
determinants of Toeplitz matrices 0, S [HI as 



G-i G-2 

Go G_i 

Gr-2 Gr-3 

Gi Go 
G2 Gi 

Gj' Gi'—i 
{a'-f - GrG. 



• ■ G-r^l 

■■ G_i 

G-r+2 
G-r+3 

Gi 



(18) 



where 

Gr 



)tanh ( ^2'^ ] cos{r4 



[1 + Acos( 
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/ d(j) tanh 







sm{rq 



UJ4, 

sin (h 
) 

W0 



(19) 



with r = |i — j| being the distance between the two spins 
at the sites i and j (for the nearest- neighbour (n.n.) case, 
r = l). 

We next consider the interaction of the chain of qubits, 
each qubit representing a S ~ ^ spin in transverse-field 
XY model, with an environment. We choose the initial 
state (time t = 0) oi the whole system at time t = to 
be of the product form, i.e., 

p(0) =/9,(0)®Pe(0) (20) 

where the density matrices ps and pe correspond to the 
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Figure 3. (Color online) (a) Variations of |ci(0)| and |c2(0)| 
with A for 7 — 0.7 (|ci(0)| solid line and |c2(0)| dashed line) 
and 7 — 0.8 (|ci(0)| dotted line and |c2(0)| dot-dashed line), 
(inset) Variation of (= |c2(0)| / |ci(0)|) with A for 7 = 0.7 
(solid line) and 7 = 0.8 (dot-dashed line), (b) Variations 
of |ci(0)| and |c2(0)| with 7 for A = 0.7 (|ci(0)| dotted 
line and |c2(0)| dot-dashed line) and A = 0.8 (|ci(0)| solid 
line and |c2(0)| dashed line), (inset) Variation of Rc with 
7 for A — 0.7 (dot-dashed line) and A — 0.8 (solid line). 
-Rc = |c2(0)| / |ci(0)| for 7 > (BPF channel) and Rc = 
|ci(0)! / |C2(0)| for 7 < (BF channel). 



system and environment respectively. We assume that 
the environment is represented by L independent reser- 
voirs each of which interacts locally with a qubit consti- 
tuting the XY chain. The reduced density matrix ob- 
tained from Eq. (|20l) can be written as 



Pr{0) = prs{0) <E) Pre{0) 



(21) 



where prs and p^e represent the two-qubit reduced den- 
sity matrix of the transverse-field XY chain and the cor- 
responding reduced density matrix of the two-reservoir 
environment respectively. The quantum channel describ- 
ing the interaction between a qubit and its environment 
can be of various types [H, . We first assume the two 
qubits to represent nearest-neighbour spins in the XY 
chain. In this paper, we investigate the dynamics of the 
two-qubit classical and quantum correlations (in the form 
of the QD) under the influence of the bit-flip (BF), bit- 
phasc-flip (BPF) and phase-flip (PF) channels. In the 
Kraus operator representation, an initial state, prs(O), of 
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the qubits evolves as HO, Ull 

Prs{t)^Y.^>^'''P-^^^)^U (22) 

where the Kraus operators E^,ij = ® Ei, satisfy the 
completeness relation i/ Ef^.iyEj^ v = I for t- The 
BF, BPF and PF channels destroy the information con- 
tained in the phase relations without involving an ex- 
change of energy. The Kraus operators for these channels 
are 

^0 ^ ^ ( 1 ) ' ^1 = ^'^'^ (2^) 

where i = a; for the BF. « = y for the BPF and i = z 
for the PF channel with q' ~ 1 — p/2. In the case of 
Markovian time evolution, p is given by p = 1 — e~^* 
with 9 denoting the decay rate. A fuller description of 
the channels and the corresponding dynamics can be ob- 
tained from Refs. [13, [13, Hal- Using the two-spin reduced 
density matrix described by Eqs. © and ([TS|) - ([TO)) as the 
initial state prs(O), one can calculate the time-evolved 
state Prs{t) for a specific quantum channel using the Eqs. 

and Since Prs{t) has the same form as Eq. 

([7]), the mutual information, I{pab), the QD Q{pab) 
and the classical correlations C{pab) can be computed 
at any time t with the help of the formulae in Eqs. (|S])- 
(fn|) . The main results of our calculations for the various 
quantum channels are described in the following section. 



III. RESULTS 

A. Nearest-neighbour quantum correlations at zero 
temperature 

1. Bit-flip and bit-phase-flip channels 
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Figure 4. (Color online) (a) Variation of Pac as a function 
of A with 7 = 0.7 (solid line) and 7 = 0.8 (dot-dashed line), 
(inset) The first derivative of psc w.r.t. A diverges as the QCP 
Ac = 1 is approached, (b) Variation of psc as a function of 
7 with A — 0.7 (solid line) and A — 0.8 (dot-dashed line) for 
the BPF channel, (inset) (i) The first derivative of psc w.r.t. 
7 exists when 7c > for the BPF channel, (ii) Variation of 
Psc as functions of 7 with A = 0.7 (solid line) and A — 0.8 
(dot-dashed line) for the BF channel. 



We hrst consider the case of T = so that the XY 
spin chain is in its ground state. In the case of the BPF 
channel, the dynamical evolutions of the mutual informa- 
tion I{pab), the classical correlations C (pab), the QD 
Q (pab) and the absolute values of the coefficients ci and 
C2 as a function of the parametrized time p {p = 1 — e~^*) 
are shown in Fig. [5] for A = 7 = 0.7. The dynamics are 
similar to type (ii) where sudden changes in the decay 
rates of the classical correlations and the QD occur at a 
parametrized time instant Psc- Inset (a) shows the varia- 
tions of |ci(p)| and |c2(p)| versus p. The crossing point is 
given by psc- Inset (b) shows the time-evolutions of the 
mutual information, the classical correlations and the QD 
for A = 0.7, 7 — —0.7 which are described by type (Hi) 
dynamics, i.e., both C (pab) and Q (pab) decay mono- 
tonically. The dynamics in the case of the BF channel 
have an interesting correspondence with the dynamics of 
the BPF channel. Type (ii) and Type (Hi) dynamics are 
obtained in the case of the BPF (BF) channel for -|-ve (- 
ve) and -ve (-l-ve) values of 7 respectively. The following 
analysis provides us with an understanding of this result. 
The rules of evolution for the coefficients ci and C2 are 
obtained from Eqs. ©, d!]), ^ and ([23]) in the cases of 
the BPF and BF channels. These are given by 



BPF Channel: 



BF Channel: 



ci(p) 


=Ci(0)(l 




C2(p) 


= C2(0) 




C3(P) 


= C3(0)(l 


-Pf 




= C4(0)(l 


-P) 


ci(p) 


=Ci(0) 




C2(p) 


=C2(0)(1 


-Pf 


C3(P) 


=C3(0)(1 


-Pf 


Ciip) 


=C4(0)(1 


-P) 



From Eqs. ^ and ([TS]), the coefficients ci(0) and C2(0) 
are identified as the two-spin correlation functions (afaj) 
and {(T^cT^j) respectively (in the present study, r ~ \i—j\ = 
1). Fig. ^a) shows the variation of the absolute values 
of the coefficients ci(0) and €2(0) with A for 7 > 0. The 
general observation is that for 7 > 0, |ci(0)| > |c2(0)| for 
A e [0, 1]. Fig. MX>) shows the plots of |ci(0)| and |c2(0)| 
versus 7 for fixed A values. One finds that when 7 is 
< 0, |ci(0)| < |C2(0)| for A e [0,1]. Since |ci(0)| > |c2(0)| 
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Figure 5. (Color online) Phase-flip channel: Decay of mutual 
information / (pab) (solid line), classical correlations C (pab) 
(dashed line) and quantum discord Q (pab) (dot-dashed line) 
as a function of the parametrized time p = 1 — e^^' for A = 
0.7, 7 = ±0.7. In this case, psc = 0.173. 



for 7 > 0, Eq. ((24| . describing the rules of evolution 
for the EPF channel, indicates that in this case there is a 
certain parametrized time instant at which |ci(p)| crosses 
|c2(p)| (= |c2(0)|). This specific value turns out to be the 
same as Psc at which a sudden change in the decay rates 
of the classical and quantum correlations occur (Fig. [5]). 
From the condition |ci(p)| = |c2(p)| and Eq. (f^. one 
obtains the following expression for psc in the case of the 
EPF channel: 

p,c = 1 - y^c (26) 

where Rc = |c2(0)| / |ci(0)|. The QD (EPF channel, 
7 > 0) is given by Q2 in Eq. ([H]) for p G [0, 1] . At p = Psc, 
(ci(p) - C2(p)) or (ci(p) -I- C2(p)) changes sign depending 
on whether ci(0) and C2(0) have the same or different 
signs respectively. Noticing that z{p) = (ci(p) + C2{p)) /4 
and fip) = (ci(p) -C2(p))/4 (Eq. (O), the magnitude 
of r which involves \z\ and |/| (Eq. ((H])) changes at 
the point p = Psc, while S (pb) and S (pab) remain un- 
changed, thus bringing about the sudden change in the 
decay rates of the quantum correlations. With 7 < 0, 
|ci(0)| is < |c2(0)| so that ci{p) never crosses C2{p) and 
one obtains type (Hi) dynamics. We next compare the 
rules of evolution for the EF channel (Eq. (|25p ') with 
those of the EPF channel (Eq. One finds that 

the rules of evolution for c^^p) and c^^p) remain identical 
while those for ci {p) and C2 (p) arc interchanged. This ex- 
plains the interchange of dynamics type mentioned ear- 
lier: for 7 > 0, since |ci(0)| is > |c2(0)|, ci{p) never 
crosses C2{p) for the EF channel so that type (in) dy- 
namics hold true whereas for 7 < 0, the crossing point is 
given by Rc = |ci(0)| / |c2(0)|. The effective interchange 
in the identities of ci(0) and 02(0) across the 7 = 
(critical XX) line is understood by remembering that 
ci(0) = (afaf^i) and 02(0) = (crfcrf^,^}. The correlation 
function (afaf^^) is interchanged with (crfcrf+i) under 
the transformation 7 — > —7 (sec Eqs. (|18p and ([TO]) '). 
The variations of Rc with A and 7 are shown in the in- 
sets of Fig. [31 

Fig. m^a) shows the plots of Psc and the first deriva- 



tive of Psc w.r.t. A, (^-ff^j (inset), as functions of A 

for the EPF channel (7 > 0). The first derivative of psc 
w.r.t. A diverges as the quantum critical point (QCP) 
Ac = 1 is approached. Fig. |H[b) shows the variation of 
Psc with the anisotropy parameter 7 for the EPF chan- 
nel. As expected, Psc 7^ for 7 > 0. Eoth |ci| and \c2\ 
individually tend to zero as 7 — > for AG [0,1] and the 
ratio Rc (= |c2(0)| / |ci(0)|) ^ 1 as 7 ^ (Fig. E]). The 
inset (i) shows the variation of the first derivative of psc 

w.r.t. 7, 1^^^^ , as a function of 7 for fixed values of A. 

The first derivative oi Psc w.r.t. 7 shows a discontinuity 
at the anisotropy transition point, 7 = 0, indicating a 
QPT. Inset (ii) shows the variation of Psc as a function 
of 7 for the EF channel {psc 7^ for 7 < 0). In a recent 
study [l^ , we have investigated the Markovian evolution 
of the classical and quantum correlations for a number 
of quantum channels in the case of the transverse-field 
Ising model, a special case (7 = 1) of the transverse-field 
XY model considered in this paper. For the EPF chan- 
nel, we have pointed out that the first derivative of the 
parametrized time instant at which the sudden changes 
in the decay rates of classical and quantum correlations 
occur w.r.t. A diverges as the Ising model critical point 
Ac = 1 is approached. We have further shown that the 
EPF (EF) channel exhibits only type (ii) (type (Hi)) dy- 
namics. In the case of the XY model, positive and nega- 
tive values of the anisotropy parameter 7 allow for both 
type (ii) (7 < 0) and (Hi) (7 > 0) dynamics for the EF 
channel. Similarly, unlike in the case of the transverse- 
field Ising model, type (Hi) dynamics are possible for the 
EPF channel (7 < 0). 



2. Phase-flip channel 

The variations of the mutual information / (pab), clas- 
sical correlations C {pab) and QD Q{pab) versus the 
parametrized time p for (A = 0.7,7 = ±0.7) have been 
plotted in Fig. [5] in the case of the PF channel. The dy- 
namics is of type (ii) and there is a specific time interval 
during which the QD is greater in magnitude than that 
of the classical correlations. The classical correlations 
remain constant at a value / (pAs)p=i, the mutual infor- 
mation of the fully decohered state, in the parametrized 
time interval psc < p < 1. In the interval < p < Psc, 
the classical correlations decay with time. On the other 
hand, the quantum correlations, as measured by the QD, 
undergoes a sudden change in its decay rate at p = Psc 
and goes to zero in the asymptotic limit p — !• 1. One 
must note that the type of the dynamics in the case of 
the PF channel is unchanged under a change in sign of 
the anisotropy parameter 7, the reason of which can be 
easily understood by examining the evolution rules for 
the coefficients c,'s from Eq. (0), (ESI and (^3)) : 

Ci(p)=Ci(0)(l-p)2 
C2(p)=C2(0)(l-p)2 
Csip) = C3(0) 

C4(p) = C4(0) (27) 

Eoth the coefficients ci and C2 decay with time identi- 
cally in the case of the PF channel whereas the other co- 
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Figure 6. (Color online) (a) Variation of Pac as functions of 
A with 7 — 0.7 (solid line) and 7 = 0.8 (dot-dashed line) in 
the case of the PF channel, (inset) The first derivative of 
Psc w.r.t. A diverges as the QCP Ac = 1 is approached, (b) 
Variation of Psc 3'S functions of 7 with A = 0.7 (solid line) and 
A = 0.8 (dashed line) for the PF channel, (inset) The first 
derivative of psc as a function of 7 becomes discontinuous at 
the anisotropy transition point 7c = 0. 



cfficiciits remain independent of time. Since the change 
in sign of the anisotropy parameter 7 only interchanges 
the values of ci and C2, it can not affect the dynamics 
that the classical and quantum correlations undergo in 
the PF channel. In this case, the sudden change in the 
decay rates of Q (pab) and C (pab) bX p = psc is un- 
derstood by noting that in the interval < p < Pac, the 
QD is determined by Q{pab) = Q2 (Eq. (fH)) ) and for 
Psc < p < 1, Q{pab) = Qi (Eq. (US])). In the regime 
Psc < P < 1, C* {pab) is given by 



C{p 



AB 



S{pa)' 

+C?l0g2 



alog2 



d + b 



a , , b 

— T + ^log2 — T 

a + b a + b 

h 

&log2 



d+b 



(28) 



The matrix elements a, b and d are functions of only C3 
and C4 which are independent of time (Eq. ([211) )• From 
Eq. ([TT|) . S[pa)j being a function of C4 only, is also 
independent of p, thereby ensuring the constancy of the 
classical correlations C (pab)- The time instant psc can 
be determined as a solution of the equation Qi — Q2. 
Fig. El^a) shows the variation of psc as a function of the 
inverse of the strength of the transverse field, A. The 
first derivative of Psc w.r.t. A for a fixed value of the 





(a) 




(b) 



Figure 7. (Color online) (a) Variation of Pac as a function 
of the system parameters A and 7 for the BPF channel (b) 
Variation of Psc as a function of the system parameters A and 
7 for the PF channel. 



anisotropy parameter 7 diverges as the QPT point Ac = 1 
is approached (although the results for positive 7 values 
only have been shown, identical result can be obtained 
for negative 7 values as well) . Fig. EJb) demonstrates the 
variation of Psc versus 7 for fixed A values. The behaviour 
of Psc is symmetric across the anisotropy transition point 
7 = 0. On approaching 7 = 0, Psc goes to zero. The inset 
of Fig. Hl^b) shows the variation of the first derivative of 
Psc w.r.t. 7 at the QCP 7 = 0. The landscapes of psc 
versus the inverse of the transverse-field strength A and 
the anisotropy parameter 7 are shown in Fig. [T] 



B. Further-neighbour quantum correlations at zero 
temperature 

We have so far considered the spins in a pair to be n.n.s. 
The dynamics of the QD when the spins are separated by 
further neighbour distances can be investigated following 
the procedure outlined in Sec. [TTl The two qubit reduced 
density matrix /Osc(O) (Eq. (j22p ) now involves further- 
neighbour spin-spin correlation functions (Eqs. ^TE\\ and 
(fT9)) with r > 1). Fig. [8] shows the decay of the mutual 
information / (pab) (solid line), the classical correlations 
C (pab) (dashed line) and the QD Q (pab) (dot-dashed 
line) as a function of the parametrized time p = 1 — e~^* 
for A = 0.7, 7 = 0.7 in the cases of the BPF (r = 2 (a) 
and r = 3 (b)) and the PF (r = 2 (c) and r = 3 (d)) 
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(a) BPF Channel, y=0.7 
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Figure 8. (Color online) Decay of the mutual information 
I {pab) (solid line), the classical correlations C {pab) (dashed 
line) and the QD Q (pab) (dot-dashed line) as a function of 
the parametrized time p = 1 — e^"* for A = 0.7, 7 — 0.7 in 
the cases of the BPF (r = 2 (a) and r = 3 (b)) and the PF 
{r — 2 (c) and r = 3 (d)) channels. Insets of (b) and (d) show 
on a magnified scale the sudden changes in the decay rates of 
the correlations 



channels with r = 2 and 3 corresponding respectively to 
second and third neighbour distances between the spins 
in a pair. Comparing Fig. [S]with Figs. [5] and [5] one finds 
that the magnitude of Psc in the case of the PF channel 
is larger than that in the case of the BPF channel. Our 
computation shows that this feature is obtained when 
r = 1 irrespective of whether A is less or greater than the 
QCP Xc ~ 1- For r > 1, however, the magnitude oipsc in 
the case of the BPF channel is less than that in the case of 
the PF channel only for A < Ac. At A > Ac, the reverse 
trend is observed. Also, for the same type of channel, 
BPF or PF, the magnitude of Psc shifts towards larger 
values as r, the separation between the two qubits (spins) 
increases. Fig. IHlshows the divergence of %^ at the QCP 
Ac = 1 for r = 2 and 3 and the BPF (a) and the PF(b) 
channels with 7 = 0.7. Maziero et. al. [sO, HH have 
shown that in the case of the transverse-field XY chain, 
the QD for spin-pairs with r > 1 is able to signal a QPT 
whereas pairwise entanglement is non-existent for such 
distances. They have further demonstrated a noticeable 
change in the decay rate of the QD as a function of r as 
the system crosses the critical point A = 1. The decay 
rate is slower for A > 1, i.e., in the ferromagnetic phase. 
In the presence of system-environment interactions, the 
QD evolves as a function of time. We select three time 
points (three different p values) at which we compute the 
QD for r = 1,2,3 and 4. Figs. [TO] (a) and (c) show 
the results for the BPF and PF channels respectively for 
7 = 0.7 and A = 0.7, i.e., A < Ac = 1. Figs. [TU] (b) and 
(d) show the corresponding plots when A is > Ac with 
A = 1.1. We find that the decay of the QD as a function 
of r becomes considerably slower as the QCP Ac 1 is 
crossed in the case of both the BPF and PF channels. 



Figure 9. (Color online) The divergence of at the QCP 
Ac = 1 in the cases of r = 2 and r = 3 and for the BPF (a) 
and PF (b) channels with the anisotropy parameter 7 = 0.7 



C. Quantum correlations at finite temperatures 

The dynamics of the QD at finite temperatures {T ^ 0) 
and for different values of r arc investigated using the 
procedure described in Sec. |lll For infinite chain spin 
models like the XX Z chain (in the presence and absence 
of magnetic field) and the transverse-field XY chain, the 
thermal quantum discord (TQD) has been shown to de- 
tect the critical points of QPTs at finite temperatures 
(T^ . [3^ . [3^. In fact, the TQD is a better indicator of 
QPTs than thermodynamic quantities and different en- 
tanglement measures. Fig. [TT] displays our results when 
system-environment interactions, which give rise to the 
dynamics of quantum correlations, are taken into ac- 
count. The first (second) column in the figure pertains 
to the BPF (PF) channel. The finite temperature calcu- 
lations are carried out for different two-qubit distances, 
e.g., r = 1,2 and 3. An examination of Fig. [TT] shows 
that at finite temperature, the first derivative ^^^jf- no 
longer diverges at the QCP Ac = 1 but exhibits a maxi- 
mum which is shifted from the QCP. With increasing T 
and r, the shift is greater in magnitude and the maxi- 
mum becomes more broadened. The last row of Fig. 1111 
shows the estimated QCP at T 7^ and for r — 1,2,3, 



given by the A value at which the maximum of the ^ 



versus A plot occurs. In the case of the BPF channel, the 
estimated QCP is always lower than Ac = 1 for all values 
of r whereas in the case of the PF channel, the estimated 
QCP lies above the Ac = 1 value for r > 2. 



IV. CONCLUDING REMARKS 

The transverse-field XY model in Id and its derivative 
model, the transverse-field Ising model, are well-known 
spin models which exhibit QPTs. The models have been 
extensively studied in recent years from the perspective 
of quantum information theory to quantify the different 
types of quantum correlations like bipartite/multipartite 
entanglement and quantum discord (QD), present in the 
ground and thermal states of the system, as a function 
of different model parameters. One significant fact which 
emerges from these studies is that the quantum correla- 
tion measures can provide signatures of QPTs. Another 
interesting line of study involves probing the dynamics of 
different type of correlations, both classical and quantum, 
due to the inevitable interaction of a quantum system 
with its environment resulting in decoherence. In this pa- 
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Figure 10. (Color online) Decay of the QD as a function of 
r at three parametrized time points (p = 1 — e~^^) for the 
(a),(b) BPF and (c),(d) PF channels. For each channel, the 
decay slows down considerably as the QCP Ac = 1 is crossed. 
The p values considered are p = 0.05 (solid squares), p = 
0.075 (empty circle) and p = 0.1 (solid circle). Also A = 0.7 
((a),(c)), 1.1 ((b), (d)) and 7 = 0.7 
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per, we study the Markoviaii time evolution of two-qubit 
states in the Kraus operator formahsm due to the local in- 
teractions of the qubits with independent environments. 
The two-qubit state is described by the nearest-neighbour 
reduced density matrix of the transverse-field XY model 
ground state in Id. The quantum channels, represent- 
ing system-environment interactions, are of three types: 
bit-flip, bit-phase-flip and phase-flip. We investigate the 
dynamics of the mutual information I{pab), a measure 
of the total correlations in the system, the classical cor- 
relations C (pab) and the QD, Q (pab) as a function of 
the parametrized time p. An earlier study [21| identified 
three different types of dynamics for a different set of ini- 
tial states of which only the type (ii) and (in) dynamics 
are observed in the present study. A significant result 
of our study is to identify quantities, and as- 
sociated with the type (ii) dynamics which diverges and 
becomes discontinuous, respectively, as a quantum criti- 
cal point is approached. The discontinuity occurs as the 
critical point 7 = (A G [0,1]) of the anisotropy tran- 
sition is approached and the divergence is obtained as 
the quantum critical point A = 1 is approached. The 
two different types of critical point transitions belong to 
two different universality classes. We find an interest- 
ing correspondence between the dynamics of bit-fiip and 
bit-phasc-flip channels, namely, the dynamics type is in- 
terchanged across the 7 = 6 line. We provide an expla- 
nation for this feature and obtain an analytic expression 
for psc, the parametrized time point at which a sudden 
change in the decay rates of the classical and quantum 
correlations occur. The transvcrsc-field XY model has 
a phase diagram richer than that of the transverse-field 
Ising model. This is reficctcd in the varied dynamics pos- 
sible for a single quantum channel in the first case. The 
BPF and BF channels in the case of the transverse-field 
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Figure 11. (Color online) The variation of 
r — 1 (solid line), r — 2 (dotted line) and r — 3 (dashed line) 
at temperatures kT = 0.01,0.1 and 0.2. The first (second) 
column refers to the BPF (PF) channel. In (g) and (h) the 
estimated QCPs are plotted against kT for the BPF (g) and 
PF (h) channels respectively and for r — 1 (solid squares), 
r — 2 (empty circles) and r = 3 (solid circles) with 7 = 1. 



Ising model are characterized by their specific type of 
dynamics [l^l whereas in the case of the XY model, the 
dynamics type changes across the 7 = line. We have 
further extended our study to finite temperatures and to 
r (the distance of separation between the two qubits de- 
scribed by the reduced density matrix) beyond the n.n. 
distance. In the latter case for T = 0, one finds that 
the derivatives diverges at the QCP Ac = 1 as in 
the case of r = 1. We have further demonstrated that 
at a specific time point, the decay of the QD as a func- 
tion of r becomes significantly slower as the QCP Ac = 1 
is crossed. We have obtained the results for the BPF 
and the PF channels but the same feature holds true for 
other channels like the amplitude damping channel. We 
have checked that the decay rate decreases continuously 
as the QCP is approached. At finite temperatures, the 
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TQD is able to provide estimates of critical points asso- 
ciated with QPTs, this is true for both n.n. as weh as 
further-neighbour qubits though higher temperatures and 
larger values of r smear out the signatures. The general 
trends obtained arc similar to those derived in the case of 
isolated systems, i.e., when system-environment interac- 
tions are ignored. In all our computations, we have not, 



however, obtained any evidence of the type (i) dynamics 
mentioned in Sec. ID In the present study, we have con- 
sidered the time evolution of the quantum correlations to 
be Markovian in nature, an issue of significant interest 
would be to look for signatures of QPTs in the case of 
non-Markovian time evolution. 
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